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Four  different  two  dimensional  (perpendicular  to  the  ambient  magnetic  field),  plasma  fluid- 
type,  numerical  simulations  following  the  nonlinear  evolution  of  the  collisional  Rayleigh-Taylor 
instability  in  the  nighttime  equatorial  F region  ionosphere  have  been  performed.  Realistic  altitude 
dependent  ion-neutral  collision  frequencies,  recombination  rates,  and  ambient  electron  density 
profiles  were  used.  In  three  cases  (ESF  0,  1,  3)  the  electron  density  profile  was  kept  constant, 
with  a minimum  bottomside  background  electron  density  gradient  scale  length  L ^ 10  km,  but 
the  altitude  of  the  F peak  was  changed,  with  F peak  altitudes  at  340,  350,  and  430  km.  All  cases 
resulted  in  bottomside  growth  of  the  instability  (Spread  F)  with  dramatically  different  time  scales 
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for  development.  Plasma  density  depletions  were  produced  on  the  bottomside  with  rise  velocities, 
produced  by  nonlinear -polarization  Ex  H forces,  of  2.5,  12,  and  160  m see,  and  percentage  deple- 
tions of  16,  40,  and  85,  respectively.  In  one  case,  KSF  0,  the  bubble  did  not  rise  to  the  topside, 
but  in  KSF  1 and  3 topside  irregularities  were  produced  by  the  bubbles  (where  linear  theory  pre 
diets  no  irregularities.  1 In  these  three  cases  Spread  F could  be  described  from  weak  to  strong  In 
the  fourth  case  (ESF  2)  the  altitude  bf  the  F peak  was  350  km.  hut  the  minimum  I,  on  the  bottom- 
side  was  changed  to  5 km.  This  resumed  in  a bubble  rise  velocity  A-  23  m sec  and  a 60''.  depletion 
with  strong  bottomside  and  modepkt e topside  Spread  F and  a time  scale  for  development  between 
ESF  1 and  3.  Two  other  eases  ESF  0 and  0”  with  (leaks  at  330  and  300  km,  respectively  and 
bottomside  L a-  10  km  were  invest igated  via  linear  theory.  These  cases  resulted  in  extremely  weak 
bottomside  Spread  F and  no  Spread  F (entire  bottomside  linearly  stable),  respectively.  These 
Emulations  show  that,  under  appropriate  conditions,  the  collisions  I Rayleigh  Taylor  instability  must's 
linear  growth  on  the  bptfomside  of  the  F region.  This  causes  the  formulation  of  plasma  density 
depletions  (bubbles)- 'Which  rise  to  the  topside  (under  appropriate  conditions)  F region  by  polariza 
lion  Ex  D motion^  High  altitude  of  the  F peak,  small  bottomside  electron  density  gradient  scale 
lengths,  and  large  percentage  depictions  yield  large  vertical  bubble  rise  velocities,  with  the  first  two 
conditions  favoring  bottomside  lintvtr  growth  of  the  instability.  The  numerical  simulation  results 
arc  in  good  agreement  with  rocket  and  satellite  in  situ  measurements  and  radar  backscatter  measure- 
ments, including  some  of  the  recent  results  from  the  August  197 ! coordinated  ground  based  measure 
ment  campaign  conducted  by  DNA  at  Kwajalein. 
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NONLINEAR  EQUATORIAL  SPREAD  F: 
DEPENDENCE  ON  ALTITUDE  OF  THE  F PEAK 
AND  BOTTOMSIDE  BACKGROUND  ELECTRON 
DENSITY  GRADIENT  SCALE  LENGTH 


I.  INTRODUCTION 

In  the  past  few  years  much  effort  has  been  expended  to  describe  and  explain  equatorial 
Spread  F phenomena  both  from  an  experiments!  and  theoretical  viewpoint  (see  for  example: 
Balsley  et  al.,  1972;  Haerendel.  1974;  Dyson  et  al.,  1974;  Chaturvedi  and  Kaw,  1975a, b,;  1976; 
Hudson  and  Kennel.  1975;  Scannapieco  and  Ossakow,  1976;  Kelley  et  a!..  1976;  Woodman  and  La 
Hoz,  1976;  Morse  et  al.,  1977;  McClure  et  al..  1977;  Chaturvedi  and  Ossakow.  1977;  Zalesak  et 
a!.,  1977;  Szuszczewicz,  1978;  Ott,  1978;  Ossakow  and  Chaturvedi.  1978;  Ossakow  et  al.,  1978; 
Hudson,  1978;  Costa  and  Kelley,  1978a,b;  and  Kelley  and  Ott.  1978).  It  is  now  generally  believed 
that  the  collisional  Rayleigh-Taylor  instability  (or  possibly  the  E x B gradient  drift  instability; 
in  any  case  a fluid  type  gradient  instability)  initiates  equatorial  Spread  F in  the  bottomside 
evening  equatorial  F region  ionosphere  ( Balsley  et  al..  1972;  Haerendel,  1974;  Scannapieco  and 
Ossakow,  1976).  Several  recent  experimental  results  lend  credence  to  this  idea  ( Kelley  et  al., 
1976;  Woodman  and  La  Hoz.  1976;  McClure  et  al..  1977).  The  generation  of  this  instability  on 
the  bottomside  leads  to  the  formation  of  plasma  density  depletions  (bubbles).  These  bubbles 
were  predicted  by  the  nonlinear  numerical  simulation  studies  (in  the  collisional  Rayleigh- 
Taylor  regime,  denoted  CR-T)  of  Scannapieco  and  Ossakow  [1976)  and  have  been  observed  ex- 
perimentally ( Kellt  v et  al.,  1976;  Woodman  and  La  Hoz,  1976;  McClure  et  al.,  1977). 

The  initial  nonlinear  numerical  simulation  studies  ( Scannapieco  and  Ossakow,  1976) 
showed  how  a piasma  mode  (CR-T),  which  is  linearly  unstable  on  the  bottomside  only,  could, 
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by  nonlinear  polarization  E x B forces  in  the  equatorial  geometry,  produce  bubble  induced  ir- 
regularities on  the  topside,  where  the  mode  is  linearly  stable.  In  the  present  paper  we  give  a 
more  comprehensive  derivation  of  the  equations  used  by  Scannapieco  and  Ossakow  (1976)  and 
extend  the  simulations  in  the  CR-T  regime  by  changing  background  ionospheric  parameters. 
In  particular,  the  altitude  of  the  /"peak  is  changed  in  different  simulations,  as  is  the  bottomside 
background  electron  density  gradient  scale  length.  Both  of  these  changes  lead  to  important  new 
simulation  results.  However,  the  most  profound  effect  comes  from  changing  the  height  of  the  F 
peak  as  this  changes  the  effective  ion-neutral  collision  frequency. 

In  section  II  the  theory  is  presented;  section  III  contains  the  numerical  simulation  results 
and  discussion;  and  section  IV  gives  the  summary  and  conclusions. 

II.  THEORY 


In  this  section  we  present  the  two  dimensional  Or.  y;  see  Fig.  1)  set  of  plasma  fluid  equa- 
tions which  are  used  to  describe  the  nonlinear  evolution  of  the  collisionai  Rayleigh-Taylor 
(CR-T)  instability  for  equatorial  nighttime  /"region  conditions.  The  ambient  geomagnetic  field, 
B,  is  taken  to  be  constant  and  in  the  ; direction,  the  y axis  is  vertically  upward  (altitude),  gravi- 
ty, g,  is  in  the  negative  y direction,  the  x axis  points  westward,  and  g x Bis  toward  the  east. 
The  ambient  electron  density  profile  depicted  in  Fig.  1 shows  a steep  bottomside.  The  bottom- 
side  steepens  due  to  recombination  effects  and  electrodynamic  forces.  In  what  follows  we  will 
neglect  any  motions  or  variations  parallel  to  B.  Consequently,  all  spatial  derivatives  will  be  tak- 
en in  the  .v,  y plane.  The  basic  plasma  two-fluid  equations  describing  the  system  are 
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where  the  subscript  <»  denotes  species  (<>  is  electron,  i is  ion),  n is  charged  particle  number  den- 
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sity,  v is  velocity;  v R is  recombination  rate;  E is  the  electric  field;  g is  gravity;  q is  the  charge, 
v is  collision  frequency;  c is  the  speed  of  light;  and  m is  mass. 


It  should  be  noted  that  nu0  in  Eq.  (1)  is  the  equilibrium  value  of  the  density  such  that  in 
equilibrium  the  right  hand  side  of  Eq.  (1)  is  zero.  In  the  usual  form  the  right  hand  side  of  Eq. 
(1)  would  have  a source  minus  a loss.  The  loss  term  is  represented  by  ( Strobel  and 

McElroy.  1970)  and  we  have  set  the  source  term  equal  to  such  that  dnil(Jdt  — 0.  The 

recombination  rate,  »•#,  represents  a combination  of  the  rate  limiting  charge  exchange  reaction 

K. 

O + + 02  - 02+  + O 


(3) 


and  the  ion-molecule  reaction 
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0 + + N2  - NO  + + N 


such  that 


(4) 


i-R  - K,«(02)  + K2n(N2) 


(5) 

where  n( 02)  and  n(N2)  are  the  neutral  number  densities  of  02  and  N2  respectively.  In  Eq. 
(2)  the  temperature  and  neutral  wind  have  been  set  to  zero  (equivalently  we  are  in  a frame 
moving  with  the  neutral  wind). 


Equation  (2)  is  then  solved  for  the  electron  and  ion  velocities  as  follows.  The  left  hand 
side  of  Eq.  (2),  i.e.,  the  inertial  terms,  are  neglected  so  that  the  time  changes  associated  with 
these  terms  occur  over  a time  scale  longer  than  the  gyroperiod  or  collision  time  (equivalently 
the  frequency  associated  with  the  inertial  terms  is  small  compared  to  the  gyrofrequencies  and 
collision  frequencies).  For  the  electrons  we  take  “ 0.  wh*re  0„  - eB/mnc  and  for 

the  ions  i /,/ft,  « 1,  where  by  v , we  mean  ion-neutral  collisions  (••„).  The  equations  for  the 
electron  and  ion  velocities  are  then  given  by 
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where  in  (6)  we  have  neglected  the  electron  g x I motion,  because  compared  with  the  same 
ion  motion  the  ratio  is  mjmr  Use  is  made  of  quasi-neutrality  such  that  nf  — nt  - n and  we 
use 


(8) 


V J - 0,  J - <■/»(»,  - ve) 

Eq.  (8)  can  be  obtained  by  subtracting  the  electron  continuity  equation  from  the  ion  continuity 
equation  (Eq.  (I)  in  the  appropriate  forms).  The  electrostatic  approximation  is  made,  i.e., 
E — — V<f>  and  we  then  obtain 


- - -£  (V<6  x *)  ■ - 
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(10) 


Then  we  set  <<>  - 4>„  + d>],  where  </>„  is  the  zero  order  potential  and  d>!  is  the  induced  or  per- 
turbed potential  The  zero  order  solution  of  Eqs.  (9)  and  (10)  requires  that  - m, g/e. 

Note  that  this  equilibrium  value  of  the  potential  makes  —0  (see  Eq.  (7))  and  makes 
vco  ” - (g  x i/11, ).  Thus  the  equilibrium  current  is  carried  by  the  electrons  and  we  have 
J„  - - n0egi/ 11  We  then  neglect  any  other  zero  order  or  ambient  electric  fields  so  that  Eqs. 
(9)  and  (10)  become 
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If  vm  were  constant  then  ( B/cv,„ ) (gxi)  would  play  the  role  of  an  effective  electric  field  (see 
Ossakow  and  Chaturvedi,  1978)  and  — V<£|  is  the  induced  or  polarization  electric  field.  Equa- 
tions (11)  and  (12)  are  the  basic  equations  we  wish  to  solve  for  n and  <£,  in  the  x,  >•  plane  In 
these  equations  n0,  vin,  and  v R are  functions  of  y (altitude). 


We  may  obtain  the  linear  growth  rate  from  Eqs.  (11)  and  (12)  by  assuming  that 

/ \ . i(.k.x+k,y -<*l) 
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<#» i - <t> xe  ‘ ' 


(13) 


where  we  have  made  the  local  approximation  k » (dnjdy)  (1  /n„).  This  results  in  the  deter- 
minant set  of  equations 
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This  set  yields 
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Setting  u ; — wf  + iy  we  obtain 
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Equation  (17)  for  the  growth  rate,  y,  clearly  shows  that  one  can  have  growth  (positive  y)  only 
on  the  bottomsidc  of  the  F peak  where  bnjby  > 0.  Moreover,  growth  (instability)  occurs 
when  the  first  term  in  (17)  is  greater  than  the  second  Equation  (17)  clearly  shows  that  on  the 
topside  where  bnjby  < 0 y is  negative  (stable).  For  purely  horizontal  propagating  waves 
— 0)  the  above  discussion  still  applies  and  we  obtain 

1 d"<.  K 

y “T~~irT — vk 

no  (19) 

cu,  “ 0 

(20) 

and  we  can  identify  n0(bnjby ) _l  with  the  background  electron  density  gradient  scale  length 
L.  In  the  denominator  of  Eqs  (17)  and  (18)  the  second  term  is  small  compared  with  the  first  if 
kL  5*>  1 (neglecting  the  altitude  dependence)  which  is  just  the  local  approximation  Then 
Eqs  (17)  and  (18)  become 
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(22) 


Equation  (21)  shows  that  the  growth  rate  is  independent  of  | k\  and  maximizes  when  k'  - k* 
(horizontal  propagation).  Also  if  one  neglects  the  y dependence  of.  then  (22)  shows  that 
w,  — - (g/fmL)  (1  /kL)  or  - ( kl .)  -1  times  the  growth  rate. 


Hi.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  basic  equations  to  be  solved  are  (11)  and  (12).  Equation  (11)  is  put  into  dimension- 
less form  with  Q = njntt.  Equation  (11)  then  becomes 
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where  we  have  used  BnJBt  “0.  Equation  (12)  becomes 
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In  Eqs.  (23)  and  (24)  v and  vR  are  functions  of  altitude  (y).  The  ion-neutral  collision  fre- 
quency, v m,  is  given  by  ( Sirobel  and  McElroy , 1970) 

vm  -2,4  x 10  -JT  n„ 

(25) 

where  nn  is  the  neutral  density  in  cm  ~3,  and  T is  the  atmospheric  temperature  in  °K.  Equa- 
tion (25)  represents  the  collision  frequency  for  0 + in  0.  The  recombination  rate,  v R,  is  given 
by  Eq.  (5)  with  (McFarland  el  al.,  1973) 

*,  -2x10 -if'  M 

1 (26) 


1 1.2  x 10  ~'2  ^-1.  T<  750 0 K 


8 x 1°  ">4  . T > 750  ° K 


The  atmospheric  quantities  in  Eqs.  (5),  (25),  (26)  and  (27)  were  obtained  from  a Jacchia  1965 
model  neutral  atmosphere.  Both  v m and  v R used  in  the  numerical  simulations  are  depicted  in 

Fig.  2. 
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The  results  of  four  nonlinear  numerical  simulations  will  be  presented  The  simulations 
are  listed  in  Table  1 and  the  differences  in  the  ambient  zero  order  driving  conditions  are  clearly 
delineated  For  these  equatorial  spread  F (ESF)  cases  only  those  parameters  listed  in  the  table 
were  changed  All  other  conditions  were  the  same  from  simulation  to  simulation.  ESF  1 is 
taken  to  be  the  canonical  case  because  it  has  essentially  the  original  simulation  conditions  of 
Scannapieco  and  Os sakow  (1976).  In  ESF  O we  have  taken  the  n0(y)  protile  01  ESF  1 and 
moved  it  down  10  km.  i.e  , the  shape  is  the  same  but  the  nk,  at  any  altitude  changes  (see  Figs.  5 

and  11).  Similarly  in  ESF  3 the  profile  of  ESF  1 is  moved  up  80  km  (see  Figs  5 and  9).  In 

ESF  2 in  the  region  of  maximum  growth  rate  of  ESF  1 the  bottoniside  ambient  electron  densi- 
ty gradient  scale  length  has  been  made  smaller,  i.e..  the  bottoniside  profile  is  steeper,  (see  Figs 
5 and  7).  The  linear  growth  rates  for  these  four  cases  are  plotted  in  Fig  3.  For  completeness, 
linear  damping  rates  for  two  other  cases  ESF  O and  ESF  O"  are  displayed  in  Fig  4 ESF  O is 
the  same  profile  as  ESF  1 with  the  F peak  moved  down  20  km.  ESF  O"  is  also  the  same  profile 
as  ESF  l,  but  with  the  peak  moved  down  50  knv  In  arriving  at  the  curves  in  Figs  3 and  4,  a 
horizontal  perturbation  was  assumed,  i.e.,  <i,  <*  exp (/ ( A x v — toZ>]  and  use  was  made  of  Eq. 

(19)  with  »■  „ and  »■  R as  given  in  Fig.  2. 

All  four  nonlinear  numerical  simulations  were  performed  over  a two  dimensional  mesh 
corresponding  to  an  altitude  (.»■)  range  of  200  km,  i.e.,  the  highest  y point  minus  the  lowest  y 
point  was  200  km,  and  an  east-west  (v)  extent  of  8 km  (sec  Figs  5.  7,  9 and  11).  The  mesh 

was  such  that  there  were  102  points  in  the  y direction  and  42  points  in  the  v direction  and 

periodic  boundary  conditions  in  both  directions.  This  corresponds  to  Ay  —2  km  and  Av  —200 
meters  Equation  (23)  was  integrated  forward  in  time  using  fiux-correctcu  transport  algorithms 
( Boris  and  SooA,  1973),  while  Eq  (24)  was  solved  using  a Chebychev-iterative  relaxation  tech- 
nique (McDonald,  1977).  The  background  zero  order  ionospheric  electron  density.  n0,  em- 
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ployed  in  each  simulation  is  depicted  in  Figs.  5,  7,  9 and  11.  Superimposed  on  the  background 
density  at  t - 0 was  an  initial  cosinelike  perturbation  (the  same  for  each  simulation)  of  the 
form  (applied  over  the  entire  mesh), 
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(28) 


x > 3.2  km 


where  x is  in  km  and  x - 0 is  in  the  center  of  the  x mesh  (see  Figs.  5,  7,  9 and  11). 


ESF  I Case.  The  n0  profile  used  in  this  canonical  case  is  depicted  in  Fig.  5.  The  y extent 
of  the  mesh  runs  from  252  km  to  452  km  with  the  F peak  at  354  km.  The  shortest  bottomside 
background  electron  density  gradient  scale  lengths  L - n0  (dn„/6.v)  -l  are  10  km  and  essen- 
tially occur  in  the  range  252  to  298  km.  The  largest  linear  growth  rate  (see  Fig  3).  yM  — 1.01 
x 10  sec  and  occurs  at  an  altitude  of  298  km.  At  this  altitude  i’K  - 3.7  x 10  ~4  sec 

Figure  3 depicts  the  range  of  linear  instability  for  ESF  1 and  this  is  from  274  to  350  km  in  alti- 

tude. Below  274  km  L is  small,  but  recombination  prevails  (see  Fig.  2)  and  between  350  and 
354  L is  large  and  recombination  prevails.  Exactly  at  the  F peak  (354  km),  drtjdy  - 0.  i.e., 
L “ °°,  and  y — -v R (see  Eq.  (19)).  Of  course  above  354  km  (the  F peak  in  this  case),  the 
gradient  changes  sign  and  all  perturbations  are  linearly  stable.  Note  in  Fig.  3 that  for  all  cases 
we  have  plotted  only  the  regions  of  instability  (growth),  i.e.,  positive  y.  Between  274  km  and 
298  km  y goes  from  4.4  x 10  "5  sec  to  1.01  x 10  sec  ~l  and  e R goes  from  7.98  x 
10  4 to  3.73  x 10  ~4  sec  Between  298  km  and  350  km,  y goes  from  1.01  x 10  sec 

to  3.48  x 10  sec  _l  and  v R at  350  km  is  7.47  x 10  sec 
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Figure  5 exhibits  contour  plots  of  constant  n/n0  or  equivalently  n{/n0>  at  t - 2000,  4000, 
8000,  and  10,000  sec  in  the  development  of  ESF  1.  Each  is  overlayed  with  a plot  of  n0  (large 
dashed  curve)  as  a function  of  altitude.  At  t — 2000  sec  the  early  phase  of  the  growth  of  the 
collisional  Rayleigh-Taylor  instability  is  exhibited.  It  shows  the  formation  of  plasma  density 
enhancements  (+)  and  depletions  ( -),  depicted  by  small  dashed  and  solid  contours,  respec- 
tively. The  depletion  contour  represents  a 16%  depletion,  with  a maximum  depletion  of  27% 
within  this  contour;  whereas,  the  enhancement  contour  represents  a 19%  enhancement,  with  a 
maximum  enhancement  of  23%  within  this  contour.  Here  and  in  subsequent  Figures  7,  9,  and 
11  the  contour  plotting  is  such  that  the  first  (outer)  depletion  contour  n/n0  is  2 ~1/4  and  each 
succeeding  inner  contour  is  2 -l/2  times  the  previous  one.  Then  we  subtract  1.0  to  find  n,  /n0 
for  depletions.  For  example,  for  three  depletion  contours,  the  outermost  would  have  «|  /r»0  — 
0.84  (16%  depletion),  the  next  inner  one  »)/n0  -0.59  (41%  depletion),  and  the  innermost 
fl,/rt0  - 0.42  (58%  depletion).  For  the  enhancement  contours,  the  first  outer  contour  is  2 1/4 
and  the  succeeding  inner  ones  are  2 1/2  times  the  previous  ones.  Then,  we  subtract  1.0  to  ob- 
tain /n0  for  enhancements.  For  example,  for  three  enhancement  contours,  the  outermost 
would  have  nt/n0  —1.19  (19%  enhancement),  the  next  inner  one  nl/n0  - 1.68  (68%  enhance- 
ment), and  the  innermost  «, /n0  - 2.38  (138%  enhancement).  With  this  scheme,  for  the  de- 
pletions (enhancements)  every  other  contour  of  rt|/n0  is  a factor  of  two  smaller  (larger). 

At  t- 4.000  sec.  in  Fig.  5.  a clear  bubble  (density  depletion)  in  the  center  of  the  x direc- 
tion is  beginning  to  form  and  rise  toward  the  F peak  with  an  innetmost  depletion  contour  of 
41%  with  a maximum  depletion  within  this  contour  of  54%.  The  innermost  enhancement  con- 
tour is  68%.  with  a maximum  enhancement  of  84%  within  this  contour.  We  also  note  the  for- 
mation of  other  depictions  in  the  wings  (near  |.v|~-  4 km).  At  t -8000  sec,  the  density  deple- 
tion or  bubble  upper  boundary  is  at  the  F peak.  The  innermost  contour  of  the  rising  bubble  is 
a 41%  depletion;  however,  there  is  a 58%  depletion  contour  below  this  at  an  altitude  ~ 300  km 
and  within  this  contour  a maximum  depletion  of  66%.  In  the  wings  there  are  depletions  with 


innermost  contours  of  71%  depletion.  There  arc  plasma  density  enhancements  at  y ^ 302  km 


with  innermost  contours  of  138%  enhancements  and  maximum  enhancements  inside  these  in- 
nermost contours  of  226%. 

At  t — 10,000  sec  (as  far  as  the  ESI-'  1 simulation  was  carried),  the  main  bubble  is  clearly 
through  the  F peak  with  an  innermost  depletion  contour  of  41%  However,  in  the  ionosphere 
below  the  bubble  near  x — 0,  there  is  a 71%  depletion  contour  with  a maximum  depletion  in- 
side this  contour  of  73%.  In  the  wings  there  are  depletions  with  innermost  contours  of  71% 
and  inside  of  these  depletion  contours  maximum  depletions  of  75%  The  innermost  enhance- 
ment contour  is  236%  with  a maximum  inside  this  contour  of  244%. 


Although  there  is  a main  bubble  which  is  beyond  the  F peak  (top  of  bubble  at  375  km) 
at  this  time,  it  has  not  completely  broken  from  the  lower  altitude  region  and  is  forming  a trail- 
ing tail.  The  top  of  the  main  bubble  in  the  central  x region  is  at  an  altitude  — 375  km  with 
the  bottom  of  the  trail  at  an  altitude  ~ 270  km.  The  widest  part  of  the  top  outermost  contour 
(16%  depletion)  of  the  bubble  is  — 3 km,  whereas;  the  innermost  contour  (41%  depletion)  is 
— 0.5  km  wide.  Also  from  t — 8,000  to  t — 10 4 sec  in  Fig  5,  the  top  most  part  of  the  bubble 
has  moved  — 24  km  which  corresponds  to  a rise  velocity  — 12  nr/sec.  At  this  time  the 
enhanced  regions  encompass  a large  region  in  the  east-west  (x)  direction  although  confined  in 
altitudes  < 320  km  (basically  isolated  depletions  move  upward,  while  isolated  enhancements 
move  downward  in  the  equatorial  spread  F geometry,  see  Ossakow  ami  Chaturwdi , 1478) 


Figure  6 displays  contours  of  constant  induced  polentii  ‘ <A,,  from  Eq  (24),  for  ESF  I 
over  the  computational  mesh  at  t ~ 10 4 sec.  The  plus  (inside  solid  lines)  denotes  positive 
values  and  the  minus  (inside  dashed  lines)  denotes  negative  values  of  potential.  Zero  potential 
falls  between  the  two  in  such  a way  that  the  outermost  solid  contour  represents  a potential  of 
5 x 10  4 statvolts  (0.15  volts)  and  the  outermost  dashed  contour  represents  a potential  of 
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—5  x 10  “4  statvolts.  Each  inner  solid  contour  represents  an  increment  of  5 x 10  ~4  stat- 
volts  and  similarly  each  inner  dashed  contour  represents  an  increment  of  — 5 x 10  ~4  stat- 
volts. Consequently,  the  potential  difference  between  any  two  solid  (dashed)  contours  is 
5 x 10  — 4 ( — 5 x 10  ~4)  statvolts.  However,  the  potential  difference  between  the  outermost 
solid  and  dashed  contours  is  10  ~3  statvolts  (0.3  volts).  This  contour  level  scheme  is  the  same 
in  Figs.  8 and  10.  We  see  that  the  more  isolated  high  altitude  part  of  the  central  bubble  is  act- 
ed on  by  an  induced  electric  field  which  points  from  west  (positive  x)  to  east  (negative  x)  and 
is  dipolar  in  nature.  This  causes  the  bubble  to  rise  with  a velocity  ( — c/fl)V  <f>,  x £.  Howev- 
er, the  lower  portion  of  the  mesh  is  acted  on  by  an  induced  electric  field  which  points  from 
east  to  west  and  is  much  weaker  than  the  induced  electric  field  acting  upon  the  isolated  portion 
of  the  central  bubble.  This  basically  causes  the  enhancements  to  move  downward  along  with 
the  lower  altitude  depletions.  Thus  the  lower  part  of  the  central  bubble  becomes  "captured"  by 
the  enhancements.  As  stated  previously  isolated  depletions  and  enhancements  should  move 
up  and  down  respectively,  at  the  equator.  However,  this  concept  is  altered  when  depletions  are 
surrounded  by  enhancements  or  vice  versa. 

ESF  2 Case.  The  n0  profile  used  in  this  case  is  displayed  in  Fig.  7.  The  y extent  is  the 
same  as  ESF  1 with  the  F peak  at  354  km.  However,  in  this  case  the  bottomside  is  steeper 
than  ESF  1 and  the  shortest  bottomside  background  electron  density  gradient  scale  lengths  L 
are  — 4.8  km  and  occur  in  the  altitude  range  between  296  and  304  km.  The  largest  linear 
growth  rate  (see  Fig.  3),  yM  - 3.2  x 10  ~3  sec  -1  (a  factor  of  three  larger  than  ESF  1)  and  oc- 
curs at  an  altitude  of  304  km.  At  this  altitude  v R -3.09  x 10  H sec  The  range  of  linear 
instability  for  ESF  2 is  depicted  in  Fig.  3.  This  goes  from  282  km  to  352  km  in  altitude.  At 
the  lowest  altitudes  n0  is  approximately  the  same  for  ESF  1 and  2.  This  makes  L larger  at 
lower  altitudes  for  ESF  2.  Consequently,  while  ESF  1 is  linearly  unstable  down  to  an  altitude  of 
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274  km,  ESF  2 is  unstable  only  down  to  282  km.  Between  282  and  304  km  y goes  from  4.5  x 
10  _s  sec  -1  to  3.2  x 10  “3  sec  _l  and  v R goes  from  6.19  x 10  sec  -1  to  3.09  x 10  ~4 
sec  Between  306  km  and  352  km,  y goes  from  2 x 10  ~3  sec  -1  to  6.82  x 10  _5  sec  -l 
and  v R at  352  km  is  7.02  x 10  ~5  sec 

Figure  7 shows  contour  plots  of  constant  n/n0  (or  n,/«0)  at  t - 1000,  3000,  4000,  and 
5000  sec  in  the  nonlinear  development  of  ESF  2.  Also  depicted  in  each  snapshot  in  this  figure 
is  a plot  of  n0  as  a function  of  altitude  for  this  simulation.  At  / — 1000  sec  the  innermost  de- 
pletion contour  represents  a 58%  depletion  with  a maximum  depletion  inside  this  contour  of 
60%.  The  enhancement  contour  represents  a 19%  enhancement.  However,  because  of  the 
enhancement  contour  spacings,  there  exists  a 63%  maximum  enhancement  inside  the 
enhancement  contour.  This  first  frame  clearly  shows  that  ESF  2 is  developing  faster  than 
ESF  1.  At  t — 3000  sec  the  bubble  is  clearly  rising  toward  the  F peak.  In  the  upper  part  of  the 
bubble,  near  340  km,  the  innermost  depletion  contour  represents  a 71%  depletion.  However,  a 
maximum  depletion  of  82%  occurs  in  the  lower  portion  of  the  bubble,  near  315  km,  inside  the 
79%  depletion  contour.  At  this  time  the  enhancement  contours  are  much  more  widespread 
than  in  the  case  of  ESF  1 with  a maximum  enhancement  of  2490%.  These  enhancements  are 
confined  in  altitude  with  280  km  < y < 305  km.  At  t —4000  sec  the  main  part  of  the  bubble 
is  already  past  the  peak.  In  the  upper  part  of  the  bubble  the  innermost  depletion  contour  is 
58%;  however,  a maximum  depletion  of  68%  occurs  inside  the  third  (innermost)  depletion  con- 
tour at  lower  (~  320  km)  altitudes,  i.e.,  within  the  trail.  Some  well  defined  depletions  (bub- 
bles) are  also  evident  in  the  lower  regions  (—  310  km  altitude)  to  the  east  and  west  of  the 
main  bubble  (|x|  — - 2 km).  Also  at  this  time  there  are  large  regions  of  enhancements 
confined  to  altitudes  275  km  ^ y ^ 310  km  with  maximum  enhancements  of  2936%. 

At  t — 5000  sec  (as  far  as  ESF  2 was  run),  the  bubble  clearly  is  well  past  the  F peak  leav- 
ing a trail  in  the  lower  ionosphere.  However,  the  trail  is  not  as  well  captured  as  in  ESF  1.  The 
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bubble  and  trail  here  extend  for  about  100  km  in  altitude  as  was  the  case  in  ESF  1.  The  inner- 
most contour  of  the  high  altitude  part  of  the  bubble  (top  of  bubble  at  387  km  altitude)  is  now 
only  41%  and  there  is  a maximum  depletion  of  73%  in  the  lower  portion  of  the  bubble,  near 
325  km,  within  the  71%  depletion  contour  The  widest  part  of  the  16%  depletion  contour  near 
the  top  of  the  bubble  is  — 3 km;  whereas,  the  41%  depletion  (inner)  contour  is  — 0.5  km 
wide.  Between  t - 3000  and  4000  sec  the  top  of  the  bubble  has  moved  up  ~ 23  km  which 
corresponds  to  a rise  velocity  — 23  m/sec.  There  is  a slight  rise  in  the  bubbles  to  the  east  and 
west  of  the  main  bubble  (compare  with  t - 4000  sec).  The  enhancements  are  confined 
between  275  km  < y < 310  km  with  a maximum  enhancement  of  3130%. 

Figure  8 depicts  contours  of  constant  induced  potential  <b\  for  ESF  2 over  the  computa- 
tional mesh  at  i — 5000  sec.  The  contours  have  the  same  scheme  as  for  ESF  1.  Once  again 
this  figure  shows  that  the  more  isolated  high  altitude  part  of  the  bubble  is  acted  upon  by  an  in- 
duced polarization  electric  field  which  points  from  west  to  east  and  is  somewhat  dipolar  in  na- 
ture. This  field  is  stronger  than  in  the  case  of  ESF  1 (see  Fig.  6;  note  that  in  ESF  2 there  are 
more  potential  contours  and  the  spacing  between  the  outermost  negative  contour  and  outer- 
most positive  contour  is  closer  for  ESF  2 than  ESF  1)  and  accounts  for  the  more  rapid  rise 
(through  —(c/B)  V<f>|  x £).  The  lower  portion  of  the  mesh,  like  ESF  1,  is  acted  upon  by  an 
induced  electric  field  which  points  from  east  to  west  and  is  much  weaker  than  the  induced 
electric  field  which  acts  on  the  isolated  portion  of  the  central  bubble.  This  field  keeps  the 
lower  portion  containing  enhancements  and  depletions  si  iow  altitudes. 

ESF  3 Case.  Displayed  in  Fig  9 is  the  nlt  profile  used  in  this  case  The  v extent  is  now 
between  332  km  and  532  km  with  the  F peak  at  434  km  This  F peak  is  80  km  higher  than 
ESF  1 and  is  about  at  the  altitude  that  Kelley  ei  a I |1976|  measured  the  F peak  after  en- 
countering a bubble  below  the  peak.  The  basic  shape  of  the  profile  is  the  same  as  F.SF  I so 
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t hut  ihe  shortest  boltomside  background  electron  density  gradient  scale  lengths  L — 10  km 
and  occur  in  the  altitude  range  between  332  km  and  378  km  (note  this  is  just  80  km  higher 
than  ESF  I)  The  largest  linear  growth  rate  (see  Fig.  3),  y —6.07  x 10  sec  -1  and  occurs 
at  an  altitude  of  378  km  where  i ■ K — 3.2  x 10  _s  sec  This  linear  growth  rate  is  a factor  of 
six  larger  than  the  maximum  growth  rate  in  ESF  1.  Unlike  ESF  1 or  2,  in  ESF  3 the  entire 
boltomside  is  linearly  unstable.  Between  332  km  and  432  km  the  recombination  term  in  Eq. 
(19)  does  not  balance  the  gradient  term.  In  this  altitude  range  r K goes  from  1.3  x 10  4 
sec  -l  to  6.4  x 10  ~6  sec  -1  while  y goes  from  2.5  x 10  sec  to  1 x 10  4 sec  We 
also  note  again  at  the  F peak,  434  km,  L — 00  and  there  is  damping  of  the  instability. 


Figure  9 displays  contour  plots  of  constant  n/n0  (or  n\/n„)  at  t — 300,  700,  1000,  and  1400 
sec  in  the  nonlinear  evolution  of  ESF  3.  Once  again  n„  as  a function  of  altitude  for  this  simu- 
lation is  depicted  in  each  of  the  frames.  At  r — 300  sec  the  early  phase  is  exhibited  with  a 16% 
depletion  contour  with  a maximum  depletion  inside  this  contour  of  26%.  The  enhancement 
contour  represents  a 19%  enhancement  with  a maximum  23%  enhancement  inside  this  con- 
tour. Clearly  this  simulation  is  developing  more  rapidly  than  either  ESF  1 or  2.  At  t — 700  sec 
a rising  bubble  in  the  central  v region  is  forming  with  an  innermost  depletion  contour  of  79%, 
with  a maximum  depletion  inside  this  contour  of  84%.  Depletion  contours  of  16%  in  the  wings 
are  also  in  evidence.  The  innermost  enhancement  contour  represents  a 138%  enhancement 
with  a 187%  enhancement  maximum  inside  this  contour.  At  t — 1000  sec  we  see  the  bubble 
rising  past  the  F region  peak.  At  this  time  the  innermost  depletion  contour,  in  the  high  alti- 
tude portion  of  the  bubble  (near  430  km),  represents  an  85%  depletion  and  there  is  a max- 
imum depletion  inside  this  contour  of  85.2%.  There  is  a long  trail  associated  with  the  bubble  at 
this  time  which  extends  about  100  km  in  altitude  However,  the  trail  is  not  as  well  captured  as 
in  FSF  1.  I he  depletions  in  the  wings,  near  | x | — 4 km,  exhibit  an  innermost  contour  of  58% 
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depletion . The  enhanced  regions  at  this  time  show  an  innermost  contour  of  572%  enhance 
ment  with  a maximum  enhancement  of  674%  inside  of  this  contour.  Here  we  also  note  that 
the  enhancements  encompass  a sizeable  region  in  the  x direction,  but  are  confined  to  altitudes 
< 390  km 

At  i — 1400  sec  (as  far  as  ESF  3 was  carried),  the  top  of  the  main  bubble  in  the  central 
x region  is  at  an  altitude  — 500  km  which  is  well  beyond  the  F peak  altitude  of  434  km. 
Moreover,  this  bubble  has  a long  trail  connecting  it  to  an  altitude  of  357  km.  The  widest  part 
of  the  top  outermost  16%  depletion  contour  is  ~ 3 km;  whereas,  the  innermost  58%  depletion 
contour  is  — 0.5  km  wide.  Therefore,  this  shows  a long  narrow  bubble  which  extends  — 150 
km  in  altitude  and  a few  km  in  the  east-west  dimension.  Radar  backscatter  observations 
( Woodman  and  La  Ho:.  1976)  have  exhibited,  on  occasion,  "plumes"  to  extend  over  100  km  in 
altitude.  They  have  identified  these  plumes  with  regions  of  depleted  plasma  density.  Within  the 
innermost  bubble  contour  at  high  altitudes,  in  Fig.  9 at  1400  sec,  there  is  a maximum  depletion 
of  70%.  The  innermost  contour  in  the  depleted  wings  represents  a 96%  depletion  with  a max- 
imum depletion  inside  this  contour  of  97%.  The  innermost  contour  in  the  enhancement  region 
corresponds  to  a 2590%  enhancement  with  a maximum  enhancement  of  2640%  inside  this  con- 
tour. Between  t — 1000  sec  and  1400  sec  the  top  part  of  the  bubble  rose  ~ 65  km  and  this 
represents  a rise  velocity  — 160  m/sec.  This  is  large  compared  with  the  12  m/sec  of  ESF  1 and 
the  23  m/sec  of  ESF  2. 

Figure  10  depicts  contours  of  constant  induced  potential,  «/> , , over  the  ESF  3 computa- 
tional mesh  at  t — 1400  sec.  Once  again  the  more  isolated  high  altitude  part  of  the  central 
bubble  (see  Fig.  9)  is  acted  on  by  a dipolar  like  induced  electric  field  pointing  from  west  to 
east.  This  causes  the  bubble  to  rise  with  the  - (c/S)V</> ] x £ velocity.  The  outermost  positive 
and  negative  potential  contours  are  more  closely  spaced  than  in  either  ESF  1 or  2 and  results 
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in  a more  rapid  rise  velocity.  The  lower  portion  of  the  mesh  is  acted  on  by  an  induced  electric 
field  which  points  from  east  to  west  and  this  is  much  weaker  than  the  higher  altitude  electric 
field  pointing  from  west  to  east.  This  causes  the  low  altitude  enhancements  and  depletions  to 
remain  at  low  altitudes  and  results  in  a capturing  of  the  central  portion  of  the  rising  bubble 

ESF  O Case.  Figure  11  depicts  the  nn  profile  used  in  this  simulation  The  y extent  is  now 
between  242  km  and  442  km  with  the  Fpeak  at  344  km  This  is  the  same  shape  profile  as  F.SF 
1 except  that  the  profile  has  been  moved  down  10  km  Consequently,  the  shortest  bottomside 
background  electron  density  gradient  scale  lengths  L — 10  km  and  occur  in  the  altitude  range 


between  242 

km  and 

288  km.  The 

largest  linear  growth 

rate 

(see 

Fig.  3), 

y \1  ™ 6 . 1 8 x 

10  “4  sec  _1 

1 and  occurs 

at  an  altitude  of 

288 

km 

where 

v n *"5.12  x 

10  “4  sec  -1 

. This  linear  growth  rate  is  about  60%  of  the  ESF  1 

maximum  linear 

growth  rate.  From  Fig.  3 we  note  that  the  ESF  O and  ESF  1 growth  rates  are  the  same  up  to 
288  km  and  both  are  stable  below  274  km.  In  Fig  3.  ESF  O shows  instability  between  274  km 
and  332  km.  This  is  over  a narrower  altitude  range  than  either  ESF  1.  2.  or  3 and  represents 
weak  instability.  In  this  altitude  range  y goes  from  4.55  * 10  sec  1 to  3.94 
x 10  sec  and  v R goes  from  7.98  x 10  ~4  sec  -l  to  1.3  x 10  4 sec 

Figure  11  depicts  contour  plots  of  constant  n/n0  (or  »|/n(l)  at  t — 5000,  10,000  and 
15,000  sec.  At  r — 5000  sec  there  is  a 16%  depletion  contour  with  a maximum  depletion  inside 
this  contour  of  33%.  There  is  a 19%  enhancement  contour  with  a maximum  enhancement  of 
33%  inside  this  contour.  This  shows  a much  slower  development  than  ESI ' 1 At  t - 10,000 
sec  a central  bubble,  disconnected  from  the  lower  depletion,  appears  to  be  rising  toward  the  7 
peak.  The  depletion  contour  of  this  bubble  is  16%  with  a maximum  depletion  of  24%  inside 
this  contour.  There  are  16%  depletion  contours  in  the  wings  at  290  km  altitude  near  | v|  — 4 
km  which  have  37%  maximum  depletions  inside  of  them  There  is  also  a wide  19%  enhance- 
ment contour  with  maximum  enhancements  inside  the  contour  of  58%. 


L»' 1 


At  t — 1 5 .OtK>  sec  the  bubble  has  risen  very  slowly  and  appears  io  be  drying  up  At  this 
time  the  top  of  the  rising  bubble  is  at  an  altitude  of  320  km  and  has  the  16%  depletion  contour 
as  at  t — 10.000  sec,  except  over  a smaller  area  There  is  a maximum  depletion  of  29%  inside 
the  lower  bubble  (near  283  km'  which  also  has  a 16%  depletion  contour  In  the  wings,  near 
290  km  altitude,  there  are  41%  depletion  contours  with  maximum  depletions  inside  these  con- 
tours of  44%.  Inside  the  19%  enhancement  contour  is  a maximum  enhancement  of  43%.  In  go- 
ing Irom  i — 10.000  sec  to  15,000  sec  both  enhancements  and  depletions  appear  to  be  decay- 
ing. even  though  the  main  bubble  is  rising  At  an  altitude  of  320  km  i #1  ~ 5000  sec  Conse- 
quently. in  this  ESF  O case,  recombination  plays  a great  role,  especially  considering  the  fact 
that  the  linear  growth  rates  are  not  very  large  Between  10.000  and  15,000  sec  the  top  of  the 
rising  bubble  has  moved  12 .5  km  This  corresponds  to  a vertical  rise  velocity  of  2.5  m/sec 
Moreover,  this  16%  depletion  contour  in  this  rising  bubble  is  — 0.7  km  wide  and  is  — 6 km 
long  in  the  altitude  direction  l or  F.SF  O,  the  enhancements  and  depletions  r.re  weak  and.  for 
the  most  part,  reside  in  an  altitude  layer  where  268  km  ^ y ^ 300  km  (bottomside)  In  this 
layer,  the  bottomside  gradient  is  sharp.  This  type  of  Spread  F condition,  i.e,  FSF  O,  could 
correspond  to  the  weak  bottomside  Spread  F identified  by  ffixx/mo/i  ami  La  Ho:  1 1 97t»l. 

LSL  O and  FSF  O"  Cases  Although  we  have  not  performed  numerical  simulations  for 
these  two  cases,  we  can  estimate  the  effects  based  on  the  four  simulations  presented  For  FSF 
O',  the  n0  profile  of  ESF  1 was  moved  down  20  km,  i.e.,  the  F peak  was  at  334  km  (10  km 
lower  than  ESF  O).  Here  the  shortest  bottomside  electron  density  gradient  scale  lengths,  L - 
10  km  and  occur  over  an  altitude  range  between  232  km  and  278  km.  For  this  case  there  is  a 
very  narrow  region  between  274  km  and  280  km  where  linear  instability  occurs  (see  Fig  41 
The  maximum  growth  rate,  y w,  is  2.15  x 10  4 sec  1 and  occurs  at  an  altitude  of  278  km 
This  is  almost  an  order  of  magnitude  smaller  growth  rate  than  FSF  I and  one-third  the  growth 
rate  of  FSF  O.  From  232  km  to  272  km  the  damping  rate  — y goes  from  2,82  x 10  'sec  1 
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to  4 26  * 1(1  ' sec  1 (see  Km  4).  From  282  km  to  4.12  km  (he  damping  rate  goes  from  2.27 
* 10  4 sec  1 to  1 27  x 10  1 sec  1 (also  see  Fig  4).  In  ESF  O (here  was  a 60  km  region 
of  linear  instability  (growth)  and  here  there  is  only  a 6 km  region  of  instability  with  smaller 
growth  rates  Consequently,  the  time  scale  for  development  here  would  he  even  longer  and 
the  strength  even  smaller  Consequently,  we  would  expect  extremely  weak  spread  F for  FSF 
O 


In  the  case  of  ESF  O"  the  «(1  profile  of  ESF  1 was  moved  down  50  km.  i.e.,  the  F peak 
was  lowered  to  304  km  Consequently,  the  shortest  L — 10  km  are  in  the  range  from  202  km 
to  248  km  However,  in  the  entire  range  from  202  km  to  402  km  altitude,  there  is  no  linear 
instability  In  this  range  the  damping  rale  (see  Fig  4)  goes  from  8.75  x 10  ’ sec  ' to  7.87 

x 10  4 sec  1 Consequently,  no  spread  F would  occur  for  this  case  ESF  O". 

IV.  SUMMARY  ANl)  CONCLUSIONS 

We  have  presented  the  results  of  four  nonlinear  numerical  simulations  of  the  collisional 
Rayleigh-Taylor  instability  for  equatorial  spread  F geometry,  using  realistic  altitude  profiles  for 
ion -neutral  collision  frequency,  recombination  rates  and  background  electron  density,  in  sec- 
tion III  Two  other  cases  where  the  altitude  of  the  F peak  was  lowered  with  respect  to  our 
canonical  case  ESF  1 were  investigated  via  linear  theory  and  conclusions  were  drawn  based  on 
the  four  numerical  simulations.  A summary  of  these  results  is  presented  in  Table  2 We  find 
that  under  favorable  conditions,  eg.  high  altitude  of  the  F peak  (small  effective  ion-neutral 
collision  frequency)  and/or  steep  bottomside  background  electron  density  gradients,  the  colli- 
sional Rayleigh-Taylor  instability  causes  linear  growth  on  the  bottomside  of  the  F region  This 
in  turn  causes  plasma  density  depletions  or  bubbles  to  be  formed  on  the  bottomside  which 
then  steepen  on  their  top  (see  Ossakow  ami  i'haiurvt'Ji,  1478)  and  nonlinearly  rise  to  the  top- 
side by  polarization  (induced)  F.  x R motion.  This  produces  irregularities  on  the  topside 


where  a linear  analysis  would  predict  no  irregularities.  High  altitude  of  the  F peak  (hp,  sec 
Table  2,  ESF  3).  small  bottomside  background  electron  density  gradient  scale  lengths  (L,  see 
Table  2,  ESF  2).  and  large  percentage  depletions  (*,/«„)  yield  large  vertical  bubble  rise  veloci- 
ties. I H.  Indeed,  we  have  shown  (sec  Table  2 and  Section  III)  that  changing  the  altitude  of 
the  F peak  from  300  km  to  430  km  can  have  dramatic  effects  on  the  evolution  of  equatorial 
Spread  F. 

With  relation  to  analytic  studies  of  rising  constant  density,  constant  shape  isolated  equa- 
torial Spread  F bubbles  in  the  collisional  Rayleigh-Taylor  regime,  (hsakow  and  Chaturvedi 
119781  have  shown  that  the  rise  velocity  is  given  by 


> « * + 


U "l 

i + 4 i — L 

b n„ 


w here  a is  one  axis  of  the  ellipse  in  the  direction  of  g,  b is  the  other  axis  in  the  east-west  direc- 
tum, and  'i[/n„  here  is  the  fractional  density  depletion  (a  positive  number  in  this  formula). 
Applying  Eq  (29)  to  the  bubble  velocities  and  n\/nt,  in  Table  2,  along  with  the  appropriate 
value  of  (Fig  2),  we  find  that  for  ESF  0,  1.  2.  and  3 1 < ti/b  < 2.  This  shows  that 

in  trying  to  model  a steepening  bubble,  which  has  a density  distribution  inside  of  it,  with  a 
constant  density-constant  shaped  bubble  results  in  the  preference  of  circular-like  shapes. 
Furthermore,  the  greater  bubble  velocity  of  ESF  3 over  ESF  1 can  be  understood  as  follows. 
Equation  (29)  and  the  general  results  of  Ossakow  and  Chaturwdi  [1978]  show  for  the  collisional 
Rayleigh-Taylor  regime 
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where  f(nx/n0)  is  an  increasing  function  of  the  fractional  depletion  nx/n0  and  depends  on 
bubble  shape  Raising  the  altitude  of  the  F peak  (as  in  ESF  3)  decreases  vm.  Also  in  the  ESF  3 
simulation  /r,  / «0  is  larger  (see  Table  2)  and  both  factors  tend  to  increase  V B.  Furthermore, 
n\ln„  is  predicted  to  be  larger  in  ESF  3 over  ESF  1,  at  least  according  to  linear  theory,  because 
nx  <x  exp  yi  where  y <*  and  is  given  by  Eq.  (19).  Raising  the  height  of  the  F peak  in- 
creases the  growth  rate,  y.  Thus,  the  ion-neutral  collision  frequency,  vin,  enters  in  two  places  in 
Eq.  (30);  whereas,  the  bottomside  background  scale  length,  L enters  only  in  / (at  least  by  the 
above  argument).  Similar  arguments  can  be  made  in  comparing  ESF  O and  2 with  ESF  1. 


To  be  sure  we  have  not  included  a horizontal  ambient  zero  order  electric  field,  E0,  in  the 
problem.  Inclusion  of  such  a field  would  in  the  E„  x B drift  frame  not  alter  Eq.  (11)  and 
would  modify  Eq.  (12)  by  adding  a term  of  the  form  E„  V (vm  n)  to  the  right  hand  side  (see 
also  Ossakow  and  Chatunedi , 1978).  This  would  also  alter  Eq.  (19),  for  example  by  changing  it 
to 


y 


, _L 

L *',*  B 


~ VR 


(31) 


where  E0  is  positive  if  in  the  eastward  direction  and  negative  if  in  the  westward  direction. 
Large  enough  values  of  a westward  £„  could,  depending  on  altitude  of  the  F peak,  result  in  sta- 
bility. Any  eastward  Elt  will  make  the  instability  grow  faster.  Similarly  an  eastward  E0  will 
make  the  bubbles  rise  faster  and  a westward  E0  will  slow  them  down  (see  Ossakow  and  Cha- 
turvedi,  1978).  In  any  case,  future  studies  including  this  effect  and  the  effect  of  a horizontal 
neutral  wind,  V (see  Ossakow  and  Chaiurvedi,  1978)  are  being  planned. 
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Our  present  studies  have  been  restrictive  in  the  east-west  (x)  dimension  to  bubble  spatial 
sizes  — 1 km  in  this  direction  (total  x dimension  — 8 km).  Our  preliminary  results  from  long 
wavelength  (■—  50-100  km)  perturbation  numerical  simulations  (Zalesak  ei  a!.,  1978)  show  that 
large  transverse  (east-west)  bubble  sizes  may  be  generated  by  the  nonlinear  evolution  of  the 
collisional  Rayleigh-Taylor  instability.  McClure  ei  at.  11977]  have  also  observed  large 
transverse  bubble  dimensions.  Future  numerical  simulation  plans  also  include  using  random 
initial  perturbations. 


Our  present  numerical  simulation  studies  as  explicitly  stated  only  covers  the  collisional 
Rayleigh-Taylor  regime.  This  allows  for  the  algebraic  solution  of  the  ion  and  electron  momen- 
tum equations  (see  Eqs.  (6)  and  (7))  and  results  in  having  to  solve  only  two  equations  (Eqs. 
(11)  and  (12))  numerically.  Including  inertial  effects  and  thus  being  able  to  handle  the  colli- 
sionless Rayleigh-Taylor  regime  requires  solving  the  full  two  dimensional  vector  ion  momen- 
tum equation  (2)  with  the  electrons  still  being  described  by  Eq.  (6).  This  is  in  addition  to  solv- 
ing a continuity  equation  (1)  and  the  divergence  of  the  current  equal  to  zero  equation  (8).  Fu- 
ture plans  call  for  making  these  changes.  However,  it  should  be  pointed  out  that  the  range  of 
validity  of  the  collisional  Rayleigh-Taylor  regime  requires  v}n  > 4j?/L  (see  for  example  Hudson 
and  Kennel , 1975;  Ossakow  and  Chaiurvedi , 1978).  Consequently,  if  we  take  L to  be  the  bubble 
scale  length,  larger  bubbles  will  remain  in  the  collisional  Rayleigh-Taylor  regime  to  higher  alti- 


The  physical  picture  of  equatorial  Spread  F that  emerges  from  our  present  simulation 
studies  can  be  described  as  follows.  After  sunset  the  E region  begins  to  recombine.  Due  to 
recombination  and  electrodynamic  effects  (or  just  recombination)  the  bottomside  F region 
background  electron  density  gradient  begins  to  steepen.  The  electrodynamic  effects  are  such 
that  an  eastward  ambient  electric  field  causes  the  bottomside  density  gradient  to  steepen  and 
the  F region  to  rise  much  akin  to  how  a barium  cloud  would  steepen  and  rise  in  equatorial 
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geometry  (see  Ossakow  and  Chaturvedi , 1978).  Thus,  the  equatorial  F region  ionosphere  can  be 
viewed  as  a giant  barium  cloud  (with  inhomogeneities  in  the  east-west  direction  occuring  over 
very  large  spatial  scale  lengths).  Indeed,  our  present  studies  of  the  nonlinear  evolution  of  the 
collisional  Rayleigh-Taylor  instability  are  similar  to  our  barium  cloud  striation  studies  (see 
Scannapteco  and  Ossakow,  1976,  Scannapieco  et  al.,  1976),  except  that  the  geometry  has  been  ro- 
tated and  gravity  substitutes  for  an  ambient  electric  field  as  the  driving  force.  This  steepening 
and  rise  can  also  be  caused  by  a downward  motion  of  the  neutral  atmosphere.  When  the  alti- 
tude is  high  enough  and/or  bottomside  background  electron  density  gradients  steep  enough  to 
overcome  recombination  effects  density  fluctuations  will  begin  to  grow  (become  linearly  un- 
stable) on  the  bottomside.  If  these  altitude  and  gradient  conditions  are  right  plasma  density 
depletions  (bubbles)  will  form  on  the  bottomside.  Similarly,  if  these  conditions  are  right  these 
bubbles  will  then  nonlinearly  rise  by  polarization  E x B motion  through  the  F peak  and  cause 
topside  Spread  F.  Although  the  present  series  of  simulations  did  not  resolve  meter  size  irregu- 
larities (in  the  present  simulations  \v  - 2 km  and  A.v  - 200  meters),  it  is  expected  that  the 
steepening  bubbles,  rising  through  the  peak,  will  bifurcate  (see  for  example  Ossakow  el  al , 
1977)  on  its  topside  and  produce  shorter  and  shorter  wavelength  irregularities  either  by  a cas- 
cade or  two  step  mechanism.  For  example,  the  gradients  associated  with  the  longer  wavelength 
<~  100  m-1  km)  unstable  fluid-type  collisional  Rayleigh-Taylor  modes  could  result  in  shorter 
wavelength  (~  1 m)  explosive  growth  kinetic  instabilities. 

Finally,  we  wish  to  compare  our  numerical  simulation  results  with  some  of  the  experi- 
mental observations.  The  phenomenon  of  rising  equatorial  Spread  F bubbles  has  been  ob- 
served by  rocket  ( Kelley  et  al.,  1976)  and  satellite  {McClure  et  al.,  1977)  in  situ  measurements 
and  inferred  from  radar  backscatter  "plume"  measurements  (Woodman  and  La  Ho: , 1976)  The 
range  of  upward  bubble  velocities,  depending  on  ambient  ionospheric  conditions,  e g height  of 
F peak  and  bottomside  gradient  scale  length,  in  our  simulations  can  account  for  many  of  the 
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experimental  observations  regarding  the  range  of  upward  velocities  reported  by  McClure  ci  al. 
(1977).  To  be  sure  a satellite  is  making  a single  point  measurement  and  does  not  have  all  the 
background  ionospheric  driving  conditions  (coordinated  measurements  with  ground  based  ex- 
periments and  possibly  rocket  in  situ  measurements  are  needed).  Nevertheless,  the  agreement 
between  the  simulations  and  the  data  (McClure  ei  al.,  1977)  is  quite  good.  The  dependence  of 
the  phenomena  on  high  ionospheres,  from  the  three  meter  backscattcr  results,  has  been  noted 
by  Farley  el  al.  (19701  and  Woodman  and  La  Ho:  [ 1 97t* j.  It  should  be  pointed  out  that  the  50 
Mllr  (3  meter  irregularities)  backscatter  is  just  a signature.  Its  full  physical  relation  to  the 
longer  wavelength  Rayleigh-Taylor  modes,  which  can  cause  Spread  F observed  on  iononograms 
and  satellite  transmission  scintillation  effects,  is  not  unequivocally  understood  at  this  point 
The  long  (in  altitude)  radar  plume  measurements  ( Woodman  and  LaHoz , 1976)  could  be  related 
to  the  long  bubbles  (e  g.,  see  F.SF  3 and  Fig.  9)  connecting  high  and  low  altitudes  (also  see  pre- 
vious discussion  in  this  section).  The  thin  layer  of  botiomside  irregularities  observed  by  W<xx/- 
man  and  La  Ho:  (19761  could  be  accounted  for  by  having  a low  ionosphere,  e.g..  our  numerical 
simulation  results  from  ESF  O and  O'.  Recently  (August  1977),  a coordinated  equatorial  F re- 
gion measurement  campaign  was  performed  under  the  auspices  of  the  Defense  Nuclear  Agency 
(DNA)  at  Kwtqalein  in  the  Marshall  Islands.  Some  of  the  preliminary  results  from  ionosonde 
(Drs.  Bibl  and  Reinisch,  University  of  Lowell),  1 meter  (150  MHz)  ALTA1R  radar  backscatter 
(Dr.  Towle,  MIT  Lincoln  Labs),  and  DNA  Wideband,  ATS-6  and  Navy  Transit  satellite  scintil- 
lation measurements  (Drs.  Baron,  Tsunoda,  Cousins  and  Livingston,  SRI)  agree  with  the 
results  of  the  present  numerical  simulations,  which  depend  on  altitude  of  the  F peak  and  steep- 
ness of  the  botiomside  background  electron  density  gradient. 
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Table  2 — Results  of  Numerical  Studies.  Note  the  bubble  rise  velocity  I #,  and  the 
fractional  depletion  njn0  refer  to  the  values  near  the  F peak  altitude 


ESF 

Altitude 

of 

F 

peak, 

/i„(km) 

Minimum 

bottomside 

Scale 

length, 

/.(km) 

Bubble 

rise 

Velocity, 

Ffl(m/sec) 

Percentage 
depletion 
n i /»„ 

Resulting 

Conditions 

1 

354 

10 

-12 

- 40% 

Maximum  linear  growth  rate 
y m — 10  3 sec  '.  Bottomside 
linear  instability  from  274  to  350 
km  altitude.  Bubble  passes 
peak  — 8000  sec.  Moderate 
top  and  bottomside  Spread  F 
with  simulation  stopped  at  104  sec. 

2 

354 

5 

-23 

-60% 

yM  — 3.2  x 10  3 sec  '.  Bottomside 
linear  instability  from  282  to 

352  km  altitude.  Bubble 
passes  peak  — 4000  sec.  Strong 
bottomside  and  moderate  topside 
Spread  F with  simulation 
stopped  at  5000  sec. 

3 

434 

10 

-160 

-85% 

y m — 6 x 10  3 sec  '.  Entire 
bottomside  from  332  to  434 
km  linearly  unstable.  Bubble 
passes  peak  — 103  sec.  Strong 
bottomside  and  topside  Spread  F 
with  simulation  stopped  at 

1400  sec. 

0 

344 

10 

-2.5 

— 1 6% 

y at  ' — 6x  10  4 sec  '.  Bottomside 
linearly  unstable  from  274  to 

332  km.  Bubble  does  not 
get  to  F peak  rather  dries 
up  due  to  recombination.  Weak 
bottomside  Spread  F only 
with  simulation  stopped  at 

15,000  sec. 

O' 

334 

10 

Undetermined 

Undetermined 

yu  — 2 x 10  4 sec  '.  Narrow 
bottomside  region  between 

274  and  280  km  linearly 
unstable.  This  would 
result  in  extremely  weak 
bottomside  Spread  F.  No 
nonlinear  simulation 
performed. 

0" 

304 

10 

None 

None 

Entire  bottomside 
between  202  and  304  km 
is  linearly  stable.  No 

Spread  F and  no  nonlinear 
simulation  performed 

z 

Fig.  I - Equatorial  Spread  F geometry.  A '(y) 
represents  the  ambient  electron  density  pro- 
flle  and  has  only  altitude  (y)  dependence. 
B points  to  the  north,  x is  west  and  so  g x B 
is  to  the  east.  A horizontal  perturbation  is 
denoted  by  k. 
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Fig.  2 - Ion-neutral  collision  frequency  (solid  line),  vJfl,  and  recombination 
rate  (dashed  line),  v^,  as  a function  of  altitude.  The  altitude  range  encom- 
passes those  altitudes  used  in  the  numerical  simulations. 


TO 


Fig.  3 - Lineal  growth  rate,  y,  versus  altitude  for  the  four  numerical  simul- 
tions  ESF  O,  ESF  I,  2.  and  3 (sec  Table  1).  A horizontal  perturbation  was 
assumed.  Only  the  regions  of  instability  are  plotted  for  these  four  cases. 


Fig.  4 - Linear  damping  rate,  y,  versus  altitude  for  ESF  O' and  O".  FSF  O 
Is  the  same  shape  as  ESF  I with  the  F peak  moved  down  20  km;  similarly 
ESF  O"  has  the  F peak  moved  down  50  km.  A horizontal  perturbation  was 
assumed.  Only  the  regions  of  stability  have  been  plotted  for  these  two  cases. 


a,  ,s‘ijUi„ViVi,i 


Itg.  5 - Contour  plots  of  constant  rt|/no  ,or  ,hc  sin'u,*,*on  ESI-  I at  t “ 2000,  4000,  8000,  anil  10,000  sec.  The  small 
tUsheil  contours  with  a plus  sign  inside  and  the  solid  contours  with  a minus  sign  inside  indicate  enhancement  and  deple- 
tions over  the  ambient  electron  number  density.  The  large  dashed  curve  depicts  the  ambient  electron  number  density 
(values  on  upper  horizontal  axis),  n^,  as  a function  of  altitude.  The  vertical  y axis  represents  altitude,  the  lower  hori- 
zontal v axis  is  east-west  range,  and  (he  ambient  magnetic  field  is  along  the  z axis,  out  of  the  figure. 


Fig.  8 - Contours  of  constant  induced  potential,  <t>  j , over  the 
ESF  2 mesh  (see  Fig.  7)  at  f = 5000  sec.  AU  other  nomen- 
clature is  same  as  Fig.  6. 
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OICY  ATTN  XRW 
OICY  ATTN  JAMES  WHELAN 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  NI CD  LIBRARY 
OICY  ATTN  ETDP  B.  BALLARD 

\ COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 

GRIFFISS  AFB,  NY  13441  |1 

OICY  ATTN  DOC  L IBRARY/TSLD 
OICY  ATTN  OCSE  V.  COYNE 

8 

L,  __  um 


SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZU 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  XPFS  MAJ  B.  STEPHAN 
01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/YA 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  YAT  CAPT  L.  BLACKWELDER 

SAMSO/SK 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS)  M.  CLAV1N 
SAMSO/MN 

NORTON  AFB,  CA  92409 
(MI NUTEMAN) 

01CY  ATTN  MNNL  LTC  KENNEDY 
COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  EEP  A.  LORENTZEN 


DEPARTMENT  OF  ENERGY 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.  0.  BOX  5400 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  FOR  A.  LABOWITZ 

EGSG,  INC. 

LOS  ALAMOS  DIVISION 

P.  0.  BOX  809 

LOS  ALAMOS,  .NM  8 5544 

0 ICY  ATTN  DOC  CON  FOR  U.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.  0.  BOX  1665 

LOS  ALAMOS,  NM  87545 

01CY  ATTN  DOC  CON  FOR  J.  WOLCOTT 

01CY  ATTN  DOC  CON  FOR  R.  F.  TASCHEK 

01CY  ATTN  DOC  CON  FOR  E.  JONES 

01CY  ATTN  DOC  CON  FOR  U.  MALIK 

01CY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

01CY  ATTN  DOC  CON  FOR  J.  ZINN 

01CY  ATTN  DOC  CON  FOR  P.  KEATON 

01CY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

SAND  I A LABORATORIES 
P.  0.  BOX  5800 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  J.  .MARTIN 

01CY  ATTN  DOC  CON  FOR  W.  BROWN 

01CY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 

01CY  ATTN  DOC  CON  FOR  T.  WRIGHT 

01CY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

01CY  ATTN  DOC  CON  FOR  3141 

01CY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 
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SAN01A  LABORATORIES 
LIVERMORE  LABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B . MURPHEY 
01CY  ATTN  OOC  CON  FOR  T.  CCXX 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 


OICY  ATTN  DOC  CON  FOR  D.  GALE 


OTHER  GOVERNMENT 


CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/S  I , RM  5G48,  HQ  BLDG 
WASHINGTON,  D.C.  20505 

0 ICY  ATTN  OSI/PSID  RM  5F  19 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20259 

CALL  CORRES:  ATTN  SEC  OFFICER  FOR) 
01CY  ATTN  R.  MOORE 

DEPARTMENT  OF  TRANSPORTATION 
OFFICE  OF  THE  SECRETARY 
TAD-49.1,  ROOM  10402-B 
400  7TH  STREET,  S.W. 

WASHINGTON,  D.C.  20590 
01CY  ATTN  R.  LEWIS 
01CY  ATTN  R.  DOHERTY 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  S INFO  ADMIN 
80ULDER , CO  80305 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01CY  ATTN  W.  UTLAUT 
0 ICY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  & ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

01CY  ATTN  R.  GRUBB 
0 ICY  ATTN  AERONOMY  LAB  G.  REID 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 

01CY  ATTN  P.  CORRIGAN 


12 


* 


m 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


r 


AEROSPACt  CORPORATION 

P.  0.  BOX  92957 

LOS  ANGELES,  CA  90009 

01CY  ATTN  I.  GARFUNKEL 
01CY  ATTN  T.  SALMI 
0 ICY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOWER 
01CY  ATTN  N.  STOCKWELL 
01CY  ATTN  D.  OLSEN 
01CY  ATTN  J.  CARTER 
0 ICY  ATTN  F . MORSE 
OICY  ATTN  SMFA  FOR  PWW 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5 OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.  0.  BOX  983 
BERKELtY,  CA  99/01 

OICY  ATTN  J.  WORKMAN 

BOEING  COMPANY,  THE 
P.  0.  BOX  5/07 
SEATTLE,  WA  98129 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 


CALIFORNIA  AT  SAN  DIEGO,  UN IV  OF 
I PAPS,  B-019 
LA  JOLLA,  CA  92093 

OICY  ATTN  HENRY  G.  BOOKER 

BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  3580/ 

OICY  ATTN  ROMEO  A.  OtLlBERlS 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

OICY  ATTN  D.  B.  COX 
OICY  ATTN  J.  P.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
b56S  ARLINGTON  BLVD 
FALLS  CHURCH,  VA  2209b 
OICY  ATTN  H.  BLANK 
OICY  ATTN  JOHN  SPOOR 
OICY  ATTN  C.  NAIL 


COMSAT  LABORATORIES 
LI  NTH  I CUM  ROAD 
CLARKSBURG,  MC  20734 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

01CY  ATTN  D.  T.  FARLEY  JR 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

01CY  ATTN  H.  LOGSTON 

01CY  ATTN  SECURITY  CPAUL  PHILLIPS) 

ESL  INC. 

*+95  JAVA  DRIVE 
SUNNYVALE,  CA  94086 

01CY  ATTN  J.  ROBERTS 
01CY  ATTN  JAMES  MARSHALL 
01CY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  S COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  ALTO,  CA  94303 

0ICr  ATTN  J.  T.  MATTINGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.  O.  BOX  8555 
PHILADELPHIA,  PA  19101 

0 ICY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COFPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

01CY  ATTN  F.  REIBERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DAS  I AC 
0 ICY  ATTN  DON  CHANDLER 
01CY  ATTN  TOM  BARRETT 
0 ICY  ATTN  TIM  STEPHANS 
01CY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 
0 ICY  ATTN  MACK  STANTON 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

HME3 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MILLMAN  14 


GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

OlCY  ATTN  JOHN  1SE  JR 
01CY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

CALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.  N.  DAVIS  (UNCL  ONLY) 

01CY  ATTN  NEAL  BROWN  CUNCL  ONLY) 

01CY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  61803 

01CY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  IL  60680 

CALL  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  J.  M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFHARD 
01CY  ATTN  JOEL  BENGSTON 

HSS,  INC. 

2 ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

INTL  TEL  8 TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAMINO  DEL  MAR 
DEL  MAR,  CA  92014 

01CY  ATTN  S.  R.  GOLDMAN 
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JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

LOCKHEED  MISSILES  S SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

01CY  ATTN  DEPT  60-12 
01CY  ATTN  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
5251  HANOVER  STREET 
PALO  ALTO,  CA  94504 

01CY  ATTN  MARTIN  WALT  DEPT  52-10 
01CY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
0 ICY  ATTN  W.  L.  IMHOF  DEPT  52-12 

KAMAN  SCIENCES  CORP 
P.  0.  BOX  7465 
COLORADO  SPRINGS,  CO  80955 
01CY  ATTN  T.  MEAGHER 

L INHABIT  CORP 
10455  ROSELLE 
SAN  DIEGO,  CA  92121 

01CY  ATTN  IRWIN  JACOBS 

LOWELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  MA  01854 

01CY  ATTN  K.  BIBL 

M.I.T.  LINCOLN  LABORATORY 
[\  P.  0.  BOX  75 

LEXINGTON,  MA  02175 

01CY  ATTN  DAVID  M.  TOWLE 
01CY  ATTN  P.  WALDRON 
01CY  ATTN  L.  LOUGHL1N 
01CY  ATTN  D.  CLARK 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.  0.  BOX  5857 
ORLANDO,  FL  52805 

01CY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5501  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.  W.  HALPRIN 
0 ICY  ATTN  TECHNICAL  LIBRARY  SERVICES 
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MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


01CY 

ATTN 

P.  FISCHER 

01CY 

ATTN 

W.  F.  CREVIER 

01CY 

ATTN 

STEVEN  L.  GUTSCHE 

01CY 

ATTN 

D.  SAPPtNF I ELD 

01CY 

ATTN 

R.  BOGUSCH 

01CY 

ATTN 

R.  HENDRICK 

01CY 

ATTN 

RALPH  KILB 

01CY 

ATTN 

DAVE  SOWLE 

01CY 

ATTN 

F.  FAJEN 

01CY 

ATTN 

M.  SCHEIBE 

01CY 

ATTN 

CONRAD  L.  LONGMIRE 

0 ICY 

ATTN 

WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.  0.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK. 

1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVE RF I ELD  BLVD. 

SANTA  MONICA,  CA  90404 

01CY  ATTN  E.  C.  FIELD  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

0 ICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  O.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.  J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

01CY  ATTN  A.  THOMPSON 
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R & 0 ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  B.  WRIGHT  JR 
01CY  ATTN  ROBERT  F.  LELEVIER 
OlCY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
OlCY  ATTN  C.  MACDONALD 
OlCY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OlCY  ATTN  CULLEN  CRAIN 
OlCY  ATTN  ED  BEDROZ IAN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OlCY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
LA  JOLLA,  CA  92038 

OlCY  ATTN  LEWIS  M.  LINSON 
OlCY  ATTN  DANIEL  A.  HAMLIN 
OlCY  ATTN  D.  SACHS 
OlCY  ATTN  E.  A.  STRAKER 
OlCY  ATTN  CURTIS  A.  SMITH 
OlCY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OlCY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 

HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

OlCY  ATTN  DALE  H.  DIVIS 

SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

OlCY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
PLEASANTON,  CA  94566 
OlCY  ATTN  SZ 


I 


SRI  INTERNAT  IONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
OICY  ATTN  L.  L.  COBB 
01CY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  & SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

OICY  ATTN  CHARLES  HUMPHREY 
OICY  ATTN  J.  W.  CARPENTER 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
UNCLASSIFIED  ONLY 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE: 

ADVANCED  RESEARCH  PROJECTS  AGENCY  (ARPA) 

STRATEGIC  TECHNOLOGY  OFFICE 
ARLINGTON,  VIRGINIA 

CAPT.  DONALD  M.  LEVINE 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

DR.  P.  MANGE 
DR.  R.  MEIER 

DR.  E.  SZUSZCZEWICZ  - CODE  7127 
DR.  TIMOTHY  COFFEY  - CODE  6700 
DR.  S.  OSSAKOW  - CODE  6780 
DR.  U.  GOODMAN  - CODE  7560 


SCIENCE  APPLICATIONS,  INC. 
1250  PROSPECT  PLAZA 
LA  JOLLA,  CALIFORNIA  92037 

DR.  D.  A.  HAMLIN 
DR.  L.  LINSON 
DR.  D.  SACHS 


DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL  LABORATORY 
NOAA 

BOULDER,  COLORADO  80302 

DR.  A.  GLENN  JEAN 
DR.  G.  W.  ADAMS 
DR.  D.  N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 


A.  F.  GEOPHYSICS  LABORATORY 
L.  G.  HANSOM  FIELD 
BEDFORD,  MASS.  01730 

DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.  M.  FORBES 
OR.  T.  J.  KENESHEA 
DR.  J.  AARONS 
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OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 

DR.  H.  MULLANEY 


COMMANDER 

NAVAL  ELECTRONICS  LABORATORY  CENTER 
SAN  DIEGO,  CALIFORNIA  92152 

DR.  M.  8LEIWEISS 
DR.  I.  ROTHMULLER 
DR.  V.  HILDEBRAND 
MR.  R.  ROSE 


U.  S.  ARMY  ABERDEEN  RESEARCH  AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MARYLAND 

DR.  J.  HEIMERL 


COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 

DR.  T.  CZU8A 


HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MASS.  02138 

DR.  M.  B.  MCELROY 
DR.  R.  LINDZEN 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PENNSYLVANIA  16802 

DR.  J.  S.  NISBET 
DR.  P.  R.  ROHRBAUGH 
DR.  D.  E.  BARAN 
DR.  L.  A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  E.  KLEVANS 


UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 

405  HILLGARD  AVENUE 

LOS  ANGELES,  CALIFORNIA  90024 

DR.  F.  V.  CORONITI 
DR.  C.  KENNEL 


UNIVERSITY  OF  CALIFORNIA,  BERKELEY 
BERKELEY,  CALIFORNIA  94720 

DR.  M.  HUDSON 


UTAH  STATE  UNIVERSITY 
4TH  N.  AND  8TH  STREETS 
LOGAN,  UTAH  84322 

DR.  P.  M.  BANKS 
DR.  R.  HARRIS 
DR.  V.  PETERSON 
DR.  R.  MEGILL 
DR.  K.  BAKER 


CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK  14850 

DR.  W.  E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 
DR.  E.  OTT 


NASA 
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